From the Kadomtsev-Petviashvili equation 
halfway to Ward's chiral model 
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Abstract 

The "pseudodual" of Ward's modified chiral model is a dispersionless limit of the matrix 
Kadomtsev-Petviashvili (KP) equation. This relation allows to carry solution techniques 
from KP over to the former model. In particular, lump solutions of the su{m) model with 
rather complex interaction patterns are reached in this way. We present a new example. 
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Ward's chiral model in 2 + 1 dimensions [Ij (see [2j for further references) is given by 

{J-^Jt)t - {J-^J.y. - {J~^Jy)y + [J'^Jx, J-^Jt] = (1) 

for an SU (m) matrix J, where Jj = dJ/dt, etc. In terms of the new variables 

xi:= {t-x)/2, X2:=y, := {t + x)/2 , (2) 

this simplifies to {J^^ Jxs)xi — {J~^ Jx2)x2 — 0) which extends to the hierarchy 

(J"'J.„+J.„-(J-V,„^J,„=0, m,n = l,2,.... (3) 

The Ward equation is completely integrabl^ and admits soliton-like solutions, often called 
"lumps". It was shown numerically [3] and later analytically [H O [6] that such lumps can 
interact in a nontrivial way, unlike usual solitons. In particular, they can scatter at right 
angles, a phenomenon sometimes referred to as "anomalous scattering"!^ Also the integrable KP 
equation, more precisely KP-I ("positive dispersion"), possesses lump solutions with anomalous 
scattering [8l[9l[10] (besides those with trivial scattering [H]). Introducing a potential (p for the 
real scalar function u via u = (px, in terms of independent variables ti,t2 (spatial coordinates) 
and (time), the (potential) KP equation is given by 

(4 - (l)titjti - 6 (j)ti - 3(7^ (j)t2t2 = , (4) 

with a = i m case of KP-I and a = 1 for KP-II. Could it be that this equation has a closer 
relation with the Ward equation? We are trying to compare an equation for a scalar with a 

^In the sense of the inverse scattering method, the existence of a hierarchy, and various other characterisations 
of complete integrability. In the following "integrable" loosely refers to any of them. 

^See also the references cited above for related work. Anomalous scattering has also been found in some related 
non-integrable systems, like sigma models, Yang-Mills-Higgs equation (monopoles) and the Abelian Higgs model 
or Ginzburg-Landau equation (vortices), see [7] for instance. 



matrix equation, and in H] the appearance of nontrivial lump interactions in the Ward model 
had been attributed to the presence of the "internal degrees of freedom" of the latter. At first 
sight this does not match at all. However, the resolution lies in the fact that the KP equation 
possesses an integrable extension to a (complex) matrix version, 

/ ti 

where we modified the product by introducing a constant N x M matrix Q, and the commutator 
is modified accordingly, so that [$t^,<^t2]Q = ^tiQ^t2 ~ ^t2Q^ti- Here $ is an M x matrix. 
If rank((5) = 1, and thus Q = VU^ with vectors U and V, then any solution of this (potential) 
matrix KP equation determines a solution (j) := U^^V of the scalar KP equation^ More 
generally, this extends to the corresponding (potential) KP hierarchies. 

Next we look for a relation between the matrix KP and the Ward equation. Indeed, there is 
a dispersionless (multiscaling) limit of the above "noncommutative" (i.e. matrix) KP equation, 

^xix3 - '^X2X2 = -f^ i^x, , ^x2]q , (6) 

obtained by introducing Xn = netn with a parameter e, and letting e — > (assuming an 
appropriate dependence of the KP variable $ on e) [2]. If rank((5) = m, and thus Q = VU^ 
with an M X m matrix U and an x m matrix V, then the m x m matrix ip := a U'^^V solves 

•Pxix^- O-"^ ^X2X2 = -[^x^-,^X2]-, (7) 

if $ solves dl]). In terms of the variables x, y, t, this become^ 

'Ptt - fxx - O"^ (fyy + [(ft - fx, fy]=0 . (8) 

Now we note that the cases a = i and a = 1 are related by exchanging x and t, hence they are 
equivalent H We choose a = 1 in the following. Then d?]) extends to the hierarchy 

<fx,^Xr,+i-fxm+ix^ = [<^x,„'PxJ, m,n=l,2,.... (9) 

The circle closes by observing that this is "pseudodual" to the hierarchy ([3]) of Ward's chiral 
model in the following sense. ([9]) is solved by 

fx„ = -J~^ Jx„+i, n = l,2,..., (10) 

and the integrability condition of the latter system is the hierarchy ([3|). Rewriting (jlOp as 
Jx„+i = ~J '■Pxni the integrability condition is the hierarchy ([9]). All this indeed connects the 
Ward model with the KP equation, but more closely with its matrix version, and not quite on 
a level which would allow a closer comparison of solutions. Note that the only nonlinearity that 
survives in the dispersionless limit is the commutator term, but this drops out in the "projection" 
to scalar KP. On the other hand, we established relations between hierarchies, which somewhat 
ties their solution structure together^ 

In the Ward model, J has values in SU{m), thus if must have values in the Lie algebra 
su{m), so has to be traceless and anti-Hermitian. Suitable conditions have to be imposed on 
$ to achieve this. Via the dispersionless limit, methods of constructing exact solutions can be 
transfered from the (matrix) KP hierarchy to the pseudodual chiral model (pdCM) hierarchy 
([9]). From [2j we recall the following result. It determines in particular various classes of (multi-) 
lump solutions of the su{m) pdCM hierarchy. 



■^See e.g. [121 113j for related ideas. 

*This Leznov equation [14' and the Ward equation arise by gauge-fixing of the hyperbolic Bogomolny equation, 
see e.g. [15] . 

^Note also that this transformation leaves the conserved density ([14[) invariant. 

^We note, however, that e.g. the singular shock wave solutions of the dispersionless limit of the scalar KdV 
equation have little in common with KdV solitons. 
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Theorem 1. Let P,T be constant N x N matrices such that = —T and = TPT ^, and 
V a constant N x m matrix. Suppose there is a constant solution K of [P,K] = -VV^T (=Q) 
such that = TRT'^. Let X be an N x N matrix solving [X,P] = 0, = TXT'^ and 
^xn+i = P"", n = 1, 2, . . .. Then Lp := —V'^T{X — K)~^V solves the su{m) pdCM hierarchy. 

Example 1. Let m = 2, N = 2, and 

-(s;.)' -(o/-)' •'K-)- 

with complex parameters a,b,c,d,p and a function / (with complex conjugate /*). Then 
^x„+i = ^xi P^ is satisfied if / is an arbitrary holomorphic function of 

w:=^x„p"-i. (12) 
n>l 

Furthermore, [P, K] = —VV'^T has a solution iff ac* + bd* = and /3 := 2S(p) ^ (where 9(p) 
denotes the imaginary part of p) . Without restriction of generality we can set the diagonal part 
of K to zero, since it can be absorbed by redefinition of / in the formula for ip. We obtain the 
following components of 

i(3 



Vii = -(/^22 = ^(|fec|2-|ad|2 + 2/39(aV)), 



9^12 = -ipli = ^[-2i{\c\^ + \d\')a*b + f3{a*df -be*!*)), (13) 

where D := (|a|2 + |6p)(|c|2 + \d\^) + |/(a;)|2 > if det(l^) 7^ 0. If / is a non-constant 
polynomial in a;, the solution is regular, rational and localized. It describes a simple lump if / 
is linear in w. Otherwise it attains a more complicated shape (see [2] for some examples). □ 

Fixing the values of 2:4, 2:5, . . ., we concentrate on the first pdCM hierarchy equation. In terms 
of the variables x,y,t given by (l2|), we then have uj = ^(t — x + 2py + p'^{t + x)), subtracting 
a constant that can be absorbed by redefinition of the function / in the solution in example 1. 
This solution becomes stationary, i.e. t-independent, \i p = ±1 The conserved density 

S := -tr[(99t - ^^f + ^y^]/2 . (14) 

of dSI) is non-negative and will be used below to display the behaviour of some solutions. 

More complicated solutions are obtained by superposition in the following sense. Given data 
{Xi, Pi,Ti,Vi) and {X2, P2,T2,V2) that determine solutions according to theorem [H we build 

The diagonal blocks of the new big matrix K will be Ki and K2. It only remains to solve 

P1K12 - K^2P2 = -ViV^T2 (16) 

for the upper off-diagonal block of K and set K21 = Tg^^ivTlgri. In particular, one can superpose 
lump solutions as given in the preceding example. 

Example 2. Superposition of two single lumps with Vi = V2 = -^2) the 2x2 unit matrix, yields 



(^11 = -^22 = -^[l32\ahi\' + (3i\ah2V + 2l3i(32'^{a*hih*2) + {(5i+(32m 



fl2 



V21 = ^[a\hi\^p2h2 + a*(3ihi\h2\^ + {b*f{a(3ihi + a*(32h2)) , (17) 
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Figure 1: Plots of 8 at times t = —90, —55, —53, 0, 30, 80 for the solution in example 2 with 
pi = — e) and p2 = i(l + e) where e = 1/20, = Aiuji, f2{uJ2) = iuj2- 




Figure 2: Origin and fate of the lump pair parts appearing in Fig. 1 (to the right in the first 
plot and to the left in the last). Plots of £: at t = -20000, -2000, 2000, 20000. 

where /3j := 29(pj), a := pi — P2, b := pi — P2, hi := a/3i/i, /12 := a*(52f2 with arbitrary 
holomorphic functions fiiyji) (where uJi is (fT2|) built with pi), respectively /2(w2), and V := 
(|6p + |/iiP)(|6p + |/i2p) + /3i/52|/ii — /i2p- This solution is again regular if pi ^ p2 ^. For 
l/il 00 (resp. I/2I 00) we recover the single lump solution ([13]) with V = I2 and / replaced 
by h (resp. /i). 

Choosing p\ = i{l — e) and p2 = i(l + e) (or correspondingly with i replaced by —i) with 
< e <C 1, and /i, /2 linear in uji, respectively uj2, one observes scattering at right angle (cf. [1] 
for the analogous case in the Ward model). 

If pi = — i(l — e) and p2 = +i(l + e), one observes the following phenomenon: two lumps 
approach one another, meet, then separate in the orthogonal direction up to some maximal 
distance, reproach, merge again, and then separate again while moving in the original direction 
[2]0 In the limit e — > 0, a vanishes and (p becomes constant (assuming /i, /2 independent of e), 
so that £ vanishes. For other choices of /i and /2 more complex phenomena occur, including a 
kind of "exchange process" described in the following. Fig. 1 shows plots of £ at successive times 
t for the above solution with /i linear in uji and /2 quadratic in uj2. The latter function then 
corresponds to a bowl-shaped lump (see the left of the plots in Fig. 2) which, at early times, 
moves to the left along the x-axis, deforming into the lump pair, shown on the right hand side of 

^See also [16] for an analogous phenomenon in case of KP-I lumps. 
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the first plot in Fig. 1, under the increasing influence of the simple lump (corresponding to the 
linear function /i) that moves to the right. When the latter meets the first partner of the lump 
pair, they merge, separate in y-direction to a maximal distance, move back toward each other 
and then continue moving as a lump pair (shown on the left hand side of the last plot in Fig. 1) 
into the negative x-direction. Meanwhile the remaining partner of the lump pair, that evolved 
from the original bowl-lump, retreats into the (positive) x-direction, with diminishing influence 
on the new lump pair, which then flnally evolves into a bowl-shaped lump (see the right of the 
plots in Fig. 2). The smaller the value of e, the larger the range of the interaction. □ 

Other classes of solutions are obtained by taking for P matrices of Jordan normal form, 
generalizing T appropriately, and building superpositions in the aforementioned sense. Some 
examples in the su{2) case have been worked out in [_2j. This includes examples exhibiting 
(asymptotic) vr/n scattering of n-lump conflgurations. The pdCM (and also the Ward model) 
thus exhibits surprisingly complex lump interaction patterns, which are comparatively well ac- 
cessible via the above theorem, though a kind of systematic classiflcation is by far out of reach. 

Acknowledgement. F M-H would like to thank the German Research Foundation for flnancial 
support to attend the workshop Algebra, Geometry, and Mathematical Physics in Goteborg. 
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